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ABSTRACT 

The 2dF Galaxy Redshift Survey has now measured in excess of 160 000 galaxy 
rcdshifts. This paper presents the power spectrum of the galaxy distribution, cal- 
culated using a direct FFT-based technique. We argue that, within the fc-space region 
0.02 < k < 0.15 /iMpc -1 , the shape of this spectrum should be close to that of the 
linear density perturbations convolved with the window function of the survey. This 
window function and its convolving effect on the power spectrum estimate are analyzed 
in detail. By convolving model spectra, we are able to fit the power-spectrum data and 
provide a measure of the matter content of the universe. Our results show that mod- 
els containing baryon oscillations are mildly preferred over featureless power spectra. 
Analysis of the data yields 68% confidence limits on the total matter density times the 
Hubble parameter f2 m /i = 0.20 ± 0.03, and the baryon fraction fi(,/Q TO = 0.15 ± 0.07, 
assuming scale- invariant primordial fluctuations. 



1 INTRODUCTION 



Present-day cosmological structure is thought to have 
formed by the gravitational amplification of small density 
perturbations. These fluctuations are readily quantified in 
terms of their Fourier modes via the power spectrum, which 
is a statistically complete description for a Gaussian field. 
The power spectrum is also of direct physical interest, be- 
cause it encodes information about the formation of the 
primordial fluctuations, and especially about how these are 
modified according to the matter content of the universe. 

In this paper, we present an estimate of the power spec- 
trum of the galaxy distribution in the 2dF Galaxy Redshift 
Survey (2dFGRS). The 2dFGRS is designed around the 2dF 



multi-fibre spectrograph on the Anglo- Australian Telescope, 
which is capable of obtaining spectra for up to 400 objects 
simultaneously over a 2 degree diameter field of view. Full 
details of the instrument a nd its performance are given in 
Lewis et al. (2000). See also http://www.aao.gov.au/2dF/. 



The survey aims to obtain redshifts for 250 000 galaxies to 
an extinction-corrected magnitude limit of 6j < 19.45. A 
description of the survey, with preliminary results, is given 
by Colless (19 99); full details of the present status c an be 
obtained from ittp : //www.mso . anu.edu.au/2dFGRS/] . 



At the time of writing, the 2dFGRS is the largest ex- 
isting galaxy redshift survey, following a natural progression 
from studies such as the CfA survey (Huchra et al. 1990), the 
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LCRS (Shectman et al. 1996), and the PSCz survey (Saun- 
ders et al. 2000). The data and analysis presented in this pa- 
per covers the sample with 166 490 redshifts observed prior 
to February 2001. A sample of this size allows large-scale 
structure statistics to be measured with very small random 
errors, and we present an initial power-spectrum analysis of 
the 2dFGRS here. Section 2 details some of the practical 
issues concerning sample selection, and Section 3 discusses 
power-spectrum estimation. The survey coverage in angular 
position and redshift is relatively complex, and the convolv- 
ing effects of the survey window are significant compared to 
the small random errors. These effects are therefore stud- 
ied in some detail, both analytically and in comparison to 
mock data, in Section 4. This leads to a robust estimate of 
the covariance matrix for the estimates of the power at dif- 
ferent wavenumbers, which is presented in Section 5. The 
covariance matrix allows proper likelihood-based model fit- 
ting, which is carried out in Section 6. The power-spectrum 
fits clearly indicate a low-density universe with Q. m h ~ 0.2, 
in agreement with many past studies. We also show that the 
preferred model requires a degree of baryon oscillations in 
the power spectrum, corresponding to a baryonic fraction 
of about 15%. We conclude by considering the consistency 
between this picture and other lines of evidence. 



2 THE 2DFGRS SAMPLE 
2.1 The angular mask 

When complete, the angular geometry of the 2dFGRS will 
consist of two declination strips plus 100 random 2-degree 
fields. One strip is near the Southern Galactic Pole (SGP) 
and covers approximately 85°xl5°; the other strip is near 
the Northern Galactic Pole (NGP) and covers 75°xl0°. 
These strips are not coplanar, which is a significant factor in 
using the survey to measure 3D structure. The 100 random 
fields are spread uniformly over the 7000 deg 2 region of the 
APM catalogue near the SGP; the present analysis includes 
71 of these fields. 

The input catalogue is a revised and extended version of 
the APM galaxy catalogue (Maddox et al. 1990a,b,c,1996). 
This includes over 5 million galaxies down to 6j = 20.5 over 
~ 10 4 deg 2 . The APM catalogue was used previously to re- 
cover the 3D power spectrum of galaxies by inverting the 
appropriate integral equations (Baugh & Efstathiou 1993; 
Efstathiou & Moody 2000). However, these techniques are 
demanding in sample variance and photometric uniformity, 
and we expect that a fully 3D analysis should yield a more 
robust result. 

An adaptive tiling algorithm is employed to cover the 
survey area with a minimum number of overlapping 2dF 
fields. With this algorithm we are able to achieve a 93% 
sampling rate with on average fewer than 5% wasted fibres 
per field. Over the whole area of the survey there are in ex- 
cess of 250 000 galaxies. At the present intermediate stage 
of the survey, unobserved fields mean that the proportion of 
targets with known redshifts is a strongly varying function of 
position. In addition, regions around bright stars are omit- 
ted, so the 2dFGRS angular mask is a complicated pattern 
on the sky (see e.g. Colless 1999). Nevertheless, because the 
tiling algorithm is known, it is possible to generate random 



catalogues that are subject to the same selection effects. 
A number of different codes have been written to achieve 
this task, with consistent results. Furthermore, because a 
3D power spectrum analysis averages over directions, small 
imperfections in reproducing the sky pattern of the real data 
tend to wash out. For example, we tried adding magnitude 
offset errors of AM — ±0.2 in each 5° Schmidt field, but 
the power spectrum did not change significantly. 

Given the sampling pattern on the sky, there are two 
possible analysis strategies: one can either build a similar 
variation into any random catalogue, or the analysis can use 
a uniform random catalogue, weighting each galaxy by the 
reciprocal of the sampling. The former strategy is superior in 
terms of shot noise, but the latter is necessary if the mask is 
correlated with real structure (e.g. fibre crowding problems 
in high-density regions). We obtain almost identical results 
with either strategy, demonstrating that the adaptive tiling 
has achieved its target of uniform selection of targets. 

2.2 Redshift selection 

The sample is chosen to be magnitude-limited at bj = 19.45 
after extinction-correcting all the magnitudes in the APM 
catalogue (Schlegel, Finkbeiner & Davis 1998). This limit 
was chosen because the mean number of galaxies per square 
degree then matches the density of fibres available with 2dF. 
The resulting distribution of galaxy redshifts n(z) dz has a 
median of approximately 0.11, and can be fitted by 

n(z) dz oc —rrdz, (1) 

[1 +(*/*)■">] 1+ ' 

where z c , 7 & 13 are fitted parameters. Fitting to all of the 
galaxy redshifts gives z c = 0.144, 7 = 2.21 & f3 = 0.554. 
However, the redshift distribution is expected to vary with 
position on the sky, because the survey depth is not com- 
pletely uniform. This arises because the spectroscopic suc- 
cess rate is a function of apparent magnitude: data from 
poorer nights are biased to brighter objects, and thus to 
lower redshifts. Also, our estimates of galactic extinction 
and CCD calibration of the zero points of the individual 
photographic plates have been revised since the original in- 
put catalogue was defined. All these effects contribute to a 
modulation of the depth of the survey, which is accounted 
for when making the random catalogue that defines the sur- 
vey volume. Because these estimates of nonuniformity can 
never be quite precise, we have chosen to allow the param- 
eters of the n(z) fit to be different in distinct zones of the 
sky, treating the NGP, SGP, and random fields separately. 
Analysis of mock catalogues shows that this makes only a 
small difference to the power estimates at k > 0.02 h Mpc -1 , 
which we use as our normal analysis limit. 



3 ESTIMATING THE POWER SPECTRUM 
3.1 The FKP method 

Following the FFT-based approach of Feldman, Kaiser & 
Peacock (1994; FKP), each galaxy in the sample was placed 
onto a 512 x 512 x 256 grid, scaled to cover the entire sample 
in each direction. The size of the cuboid that just contains all 
of the data to a redshift limit at z = 0.25 is 1388 x 1017 x 
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685 /i" 3 Mpc 3 , assuming a flat D. m = 0.3 cosmology. The 
grid used therefore leads to approximately equal Nyquist 
frequencies in each direction. For this grid, k = 0.15 h Mpc -1 
corresponds to a conservative 13% of the minimum Nyquist 
frequency. Comparing with results obtained using different 
size grids also suggests that aliasing is not a problem for 
frequencies k < 0.15, although, of course, it will affect the 
result at large k. As shown by FKP, the galaxies must be 
weighted in order to optimize the balance between cosmic 
variance and shot noise. Assuming that the fluctuations are 
Gaussian, the optimal FKP weight is 



[l + Pn(z,8, 



(2) 



where n(z, 0, <f>) is the expected galaxy density. A value of 
P = 5000 h 3 Mpc -3 was assumed, but this is not critical. 
For a high-density sample like the 2dFGRS, one very nearly 
gives equal weight to each volume element, independent of 
P. We applied an upper redshift limit of z = 0.25, to remove 
regions with a very low galaxy density where the choice of 
P would matter. We also defined a lower limit of z > 0.003. 
With restriction to objects with redshift quality flag > 3, 
this leaves a sample of 147 024 objects. 

In order to transform from redshifts to distances in 
h' 1 Mpc we need to assume values for Q m and fl„ . Strictly, 
when fitting models to the data we should alter these val- 
ues to match the model. However, the power spectrum is 
only weakly dependent on this choice: we have tried both an 
Einstein-de Sitter cosmology and a flat Q. m = 0.3 cosmology 
and find approximately the same best fit model parameters 
(see Table hi). Normally, we will present results for a flat 
flm = 0.3 cosmology. 

Before Fourier transforming the data, we need to con- 
vert from a distribution of galaxies to a distribution of over- 
densities. To do this, the survey volume was defined by a 
random catalogue that mimics the sampling of the origi- 
nal data. In order to create such a sample, we have used 
the known angular position and completeness of each field 
(described in Section 2.1), coupled with an empirically de- 
termined fit to the redshift distribution of the 2dFGRS cat- 
alogue (described in Section [^). The variations in observ- 
ing conditions, and the corresponding spectroscopic success 
rate, mean that there are in practice variations in the survey 
depth with sky position. These are allowed for by appropri- 
ate perturbations of the redshift distribution in the random 
catalogue, although this again turns out to be an unimpor- 
tant effect. The random catalogue we created had 5 times 
the number of points in the galaxy catalogue. 

Having Fourier transformed the resulting overdensity 
field, we simply subtracted the shot noise contribution from 
the result, and re-normalized P(k) in order to correct for a 
differing sample volume and FFT volume, and to correct for 
the weighting (FKP equations 2.1.3 and 2.1.10). P(k) was 
then spherically averaged over fc-space shells. 

The power spectrum of the galaxy distribution needs to 
be interpreted with care, as it is altered from that of the 
initial density fluctuations (the linear power spectrum) by a 
number of effects: 

(1) The collapse of structures. The mass no longer obeys 
the linear power spectrum on small scales after it has 
undergone non-linear collapse. 




Figure 1. The 2dFGRS window function in Fourier space. The 
fitting formula of Equation ^ is plotted (solid line) compared to 
the spherically averaged power obtained by Fourier transforming 
a random catalogue created to match the 2dFGRS window func- 
tion (solid circles). This catalogue, containing 5 times as many 
points as the true galaxy catalogue was placed ona512x512x 256 
grid that encompassed all the data. Values from individual grid 
points are also plotted (dots), and show the anisotropy of the 
window function. 



(2) Galaxies are not expected to form a Poisson sampling 
of the underlying mass fluctuations, and are biased with 
respect to this distribution. 

(3) Random oscillations of galaxies within larger collapsed 
objects causes 'fingers of God' redshift space distortions 
that damp estimates of the small-scale power. 

(4) The infall of galaxies into concentrations of mass on 
large scales (Kaiser 1987) enhances the observed power 
in the radial direction, creating a large-scale redshift 
distortion. 

Additionally, the recovered power spectrum estimated using 
the FKP approach is the convolution of the galaxy power 
spectrum with the 'window function' | Vt^fe | 2 of the survey 
(Wk being the Fourier transform of the product of the survey 
mask and redshift selection function times w(z); see FKP 
equation 2.1.7): 



PpKp(k) OC / Ptruc(k 



q)\W k (q)\ 2 d\ 



(3) 



In Section 3.2 



we show how we may model the effect of the 
window function, by convolving power spectra with an ana- 
lytic fit to the spherically averaged window function, and 
empirically by analysing realizations of Gaussian density 
fields within the survey region. The inverse problem of de- 
convolving the window function from the power spectra is 
not attempted in this work. 



3.2 The window function 

The 2dFGRS window function | Vt^fc | 2 is not compact, and 
has a significant effect on the recovered power spectrum. Al- 
though Wk is anisotropic due to the complicated real space 
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Figure 2. The effect of convolving a linear power spectrum with 
the 2dFGRS window function. Two model power spectra are com- 
pared: The top panel shows a power spectrum including a bary- 
onic component P(k,Q m h = 0.2,Qi,/Q m = 0.16) while the bot- 
tom panel shows a power spectrum P(k, Q m h = 0.2, f2(,/f2 m = 0) 
with no baryonic component. These spectra are divided by a 
smooth reference model, P(k,Q m h = 0.2,f!;,/f2 m = 0). The up- 
per solid lines show the unconvolved power spectra, and the lower 
solid lines after numerically convolving with the fit to the win- 
dow function given by Equation [l| The dashed lines show the 
average power spectrum expected from data which has had the 
average over-density artificially set to zero (see Section | S--| for 
details). Solid symbols show the shape of the average recovered 
power spectrum from 1000 Gaussian realizations of the density 
field (with different phases), placed on a 256 X 256 X 128 grid at 
locations covered by the 2dFGRS. Due to the relatively coarse 
grid used, aliasing swamps the signal at k > 0.15 h Mpc -1 , and 
data are only shown for the range of k over which the 2dFGRS 
power spectrum is fitted (Section ^j) . This region is delineated by 
the vertical dotted lines. 



selection function, we are normally interested in the spheri- 
cal average of the final power spectrum over fc-space shells. 
Therefore, if Ptruc is isotropic, we obtain the same result by 
convolving just with the spherical average of | Wk\ 2 ■ A good 
approximation to the averaged window is 

< 1 V(^ fc | 2 > = [1 + (fc/0.00342) 2 + (fc/0.00983) 4 n\ (4) 



This fit is compared to the exact 2dFGRS window in Fig. hi 
The potential cause of problems is the tail to high fe. The 
convolution involves a normalization factor J |IFfc| 2 47rfc 2 dk, 
and a significant part of this volume lies outside the 'core' 
that has a scale of 0.003 ftMpc" 1 . In this respect, the 2dF- 
GRS window is quite like an adaptive-optics PSF. 

Of course, redshift-space distortions mean that the true 
power spectrum will not be strictly isotropic in practice. 
However, the simulation results described below suggest that 
the effect of such deviations are small and do not have a 
significant impact on the recovered power spectrum. 

In order to demonstrate the effect of this window func- 
tion on power spectra, we have sampled realizations of a lin- 
ear density field covering the volume of the 2dFGRS data. 
Averaging the results of this analysis for 1000 realizations 
with different phases, and comparing with the true input 
power spectrum, we can quantify the effect of the window 
function. This is shown in Fig. g for spectra corresponding 
to two different cosmologies. The shape of the spectra are 
significantly altered, and any oscillations are damped. An- 
alytically convolving the power spectra with the fit to the 
window function given by Equation ^ provides the same re- 
sult. 

Because the random catalogue is scaled to match the 
normalization of the galaxy catalogue, the average over- 
density is artificially set to zero, forcing P(0) = 0. This 
self-normalization results in a deficit in the measured power 
spectrum equivalent to subtracting a scaled copy of the win- 
dow function, centred on k = 0. However, this is a very small 
effect in the regime of interest because |W^fc| 2 is a rapidly de- 
creasing function of k (see Fig. |l|), and P(k) is expected 
to be an increasing function of k for small k. The effect of 
this self-normalization is shown in Fig. ^| for model power 
spectra. 

3.3 Results 

The recovered 2dFGRS power spectrum is presented in 
Fig. [| To highlight features in the spectra presented in 
this paper, we have plotted their ratio with a smooth CDM 
power spectrum P(k,Q m h = 0.2, flj,/fi m = 0) that has no 
baryon features. We also take a normalization of as = 1 
for this reference model. All of the power spectra used in 
this paper are calculated using the transfer function fitting 
formulae of Eisenstein & Hu (1998) and assume a scale- 
invariant primordial spectrum unless stated otherwise. We 
have parametrized the model power spectra explicitly by 
P(k, Q m h, 0,1,/Qm) in order to avoid confusion with differ- 
ing definitions of the commonly used shape parameter F. 

The raw results are gratifyingly accurate, with frac- 
tional errors in the power of only ~ 15% out to k = 
0.02 h Mpc -1 . To within about 20%, the observed spectral 
shape is that of the P(k,£l m h = 0.2, fib/n m = 0) reference 
model between 0.02 ftMpc -1 < k < 0.6 ft Mpc -1 . At smaller 
scales, the finger-of-God redshift-space smearing is clearly 
seen to reduce the power. In many ways, the most striking 
features are the suggestions of oscillatory modulations, with 
a possible peak at k ~ 0.065 h Mpc -1 and possible troughs 
at k ~ 0.035/1 Mpc -1 and k ~ 0.1 h Mpc -1 . However, it is 
clear that the window function has caused adjacent power 
estimates to be closely correlated, so a proper covariance 
analysis is required before any significance can be given to 
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Figure 3. The 2dFGRS estimate of the redshift-space galaxy 
power spectrum, expressed as the ratio to a linear-theory CDM 
P(k,Q m h = 0.2,Q b /Q m = 0) power spectrum with n = 1 & 
eg = 1. These data do not estimate the true power spectrum, 
but give the pow er spectrum convolved with the window function 
(see Section 3.2). Error bars are determined from the diagonal 
elements of the covariance matrix (calculated in Section oh, for 
the 0.02 h Mpc" 1 < k < 0.15 h Mpc -1 data constrained by the 
vertical dotted lines. This is the region fitted in Section [i.-j. 



these apparent features. Given the precision of the basic 
power estimates, it is necessary to attain an accurate mea- 
sure of the systematic effects listed above that modify the 
shape of the spectrum. We achieve this in the next Section 
by using mock data. 



4 TESTS ON MOCK DATA 

Redshift-space and non-linear effects mean that the shape of 
the recovered power spectrum gives information about the 
linear power spectrum only at small values of k. The easiest 
way to model both these effects is via numerical simulation. 
Using an empirically motivated biasing scheme, it is possible 
to place galaxies within iV-body simulations and provide 
mock catalogues designed to mimic the 2dFGRS catalogue 
for different cosmological models (Cole et al. 1998; Baugh 
et al. 2001). In this Section we use these mock catalogues to 
demonstrate that for k < 0.15 h Mpc -1 the only important 
effect on the power spectrum is the convolution with the 
window function. 

Fig. ^| shows a collection of power spectra calculated 
from samples drawn from the ACDM and rCDM Hubble 
volu me simulations (see Cole et al. 1998, Baugh et al. 2001 



and tittp: //star-www. dur . ac . uk/~cole/mocks/main. html 



for details). The large catalogues constructed in this work 
were re-sampled depending on angular position and red- 
shift to match the correct window function for the 2dF- 
GRS data. Power spectra were calculated exactly as for 
the 2dFGRS data. Fig. [| illustrates the factors that trans- 
form the linear mass power spectrum into the non-linear 
galaxy spectrum. Panel (a) shows the power spectrum of 
the mass at z — 0. This demonstrates the increase in power 



on small scales caused by the collapse of halos. Panel (b) 
differs in that we now have to consider the effect of bias 
(artificially added to the simulation), and the window func- 
tion. In panel (c), we analyze samples designed to mimic 
the 2dFGRS data as closely as possible, including redshift- 
space effects: the finger-of-God effect that decreases small- 
scale power, and the Kaiser effect that enhances the power. 
The redshift-space and non-linear effects cancel to some ex- 
tent and give approximately the correct level of P(k) out to 
k ~ 0.5 h Mpc -1 . However, the shape of the power spectrum 
is altered for k > 0.2 h Mpc" 1 . 

We will therefore assume that, at k < 0.15/iMpc -1 , 
redshift space distortions and non-linear effects have an in- 
significant effect on the shape of P(k). Allowing the normal- 
ization to vary removes any large-scale constant bias; the 
bias is not expected theoretically to vary significantly with 
k on these large scales (Kauffmann, Nusser & Steinmetz 
1997; Benson et al. 2000). In the future, measurements of 
P ee Q.^/b from redshift -space distortions as a function of 
scale will test directly the degree to which this is true (see 
Peacock et al. 2001 for the first 2dFGRS results on redshift- 
space distortions) . This region of the power spectrum there- 
fore directly provides information about the shape of the 
linear power spectrum, and can be used with models of the 
transfer function to provide constraints on cosmological pa- 
rameters. Although this conclusion has only been justified 
here for two specific assumed models, we have performed 
similar tests on a wider range of models. In Table 1 below, 
we show explicitly how the results depend on the range of 
wavenumber considered. 



5 ESTIMATING THE COVARIANCE MATRIX 

The P(k) data points in Fig. ^| are not independent, and 
correlations extend across finite regions of fc-space. This is 
predominantly caused by the convolution with the window 
function, although there is also a contribution from non- 
linear effects (Meiksin & White 1999; Scoccimarro, Zaldar- 
riaga & Hui 1999) and redshift-space effects. Particular care 
must therefore be taken in interpreting 'wiggles' in the power 
spectrum as significant features. In order to quantify these 
correlations, we have estimated the covariance matrix for 
the data points with k < 0.15 h Mpc -1 . 

A large number of independent realizations of P(k) are 
required in order to have sufficient signal-to-noise in the co- 
variance matrix. It would be too time-consuming to perform 
separate numerical simulations for each data set. Instead, we 
have created 1000 realizations of a Gaussian random field on 
a 256 x 256 x 128 grid covering the region of the 2dFGRS 
survey. For k < 0.15 h Mpc -1 , using a smaller grid than that 
used for the 2dFGRS data does not significantly affect the 
result, and reduces the computational burden. The result- 
ing power spect ra, d etermined as for the 2dFGRS data, were 
to demonstrate the effect of convolving 



used in Section 3.2 



P(k) by the window function. These realizations provide 
an estimate of the cosmic variance within the 2dFGRS vol- 
ume. The contribution from shot noise has been calculated 
by analysing similar Monte-Carlo realizations, and has been 
included in our estimate of the covariance matrix. 

Estimating the covariance matrix in this way does not 
take into account non-linear and redshift-space effects, which 
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Figure 4. A comparison of power spectra determined from the results of the ACDM and tCDM Hubble volume simulations divided 
by a P(k,Q m h = 0.2,f2(,/C m = 0) reference model, (a) the average spectrum recovered from the z = distribution of mass in 8 cubes 
each of side 1500 h' 1 Mpc (ACDM) or 1000 h' 1 Mpc (tCDM) (solid circles with la errors: the extremal data bound the shaded area), 
compared to the scaled input power spectrum (solid line). Here, the difference between the two is dominated by the non-linear collapse 
of small-scale structures at k > 0.15 h Mpc -1 . The tCDM mass power spectrum is compared to a reference model with as = 0.6 (rather 
than as = 1.0) in order that it has approximately the same normalization as the other data, (b) the average recovered power spectrum 
from 35 (ACDM) and 30 (rCDM) real-space galaxy catalogues sampled from the simulations to match the 2dFGRS window function 
(solid circles with la errors: the extremal data bound the shaded area). These are compared to the linear spectrum convolved with 
the analytic approximation to the window function given by Equation U (solid line) . The difference between these two is dominated by 
both bias (used to create the galaxy catalogue) and non-linear effects, (c) the average recovered power spectra from 35 (ACDM) and 30 
(tCDM) redshift-space galaxy catalogues designed to fully mimic the 2dFGRS redshift-space sample (solid circles with la errors: the 
extremal data bound the shaded area), and the linear spectra convolved with the window function (solid line). Here, the differences are 



due to all of the factors listed in Section 3.1 The dashed lines i n pa nels (b) & (c) show the power spectrum expected from data in which 
the ave rage over-density is artificially set to zero (see Section for details). The 0.02 h Mpc" 1 < k < 0.15 h Mpc -1 region fitted in 



Section 6.4 is delineated by the vertical dotted lines. In calculating the error in the average power spectra measured from the simulations, 



we have assumed that the samples are independent. This is not precisely true, and the errors plotted therefore underestimate the true 
error. 
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add to the covariances. However, these effects should be 
small over the region of fc-space we are fitting. As a test 
of this we have estimated the covariance matrix using 10 
catalogues drawn from £l m h = 0.25 flb/Qm = CDM simu- 
lations with different phases by Cole et al. (1998). These 
catalogues were calculated using Cole et al. (1998) bias 
model 1. The correlations calculated over the fc-space region 
0.02/iMpc" 1 < k < 0.15 /iMpc" 1 were similar in scale to 
those calculated from the Gaussian fields. However, we find 
that the errors in P(k) calculated from numerical simula- 
tions are 16% larger than those determined from Gaussian 
simulations, although there is no evidence for a change in 
shape of the diagonal elements of the covariance matrix for 
k < 0.15 ftMpc -1 . We do see a change at k > 0.15 /iMpc" 1 
consistent with non-linear and redshift-space effects, which 
are expected to be an increasing function of k. In this work 
we adopt the conservative approach and renormalize the co- 
variance matrix calculated from the Gaussian realisations to 
match the normalisation of the numerical simulations, whilst 
keeping the correlation matrix the same. This renormaliza- 
tion does not significantly affect the primary results of this 
paper: the derived best-fit parameters are the same with or 
without this renormalization. 



6 FITTING TO THE POWER SPECTRUM 
6.1 Model parameters 

Model power spectra for different cosmologies have been cre- 
ated to compare with the 2dFGRS data using the transfer 
function fitting formulae of Eisenstein & Hu (1998). These 
formulae are essentially perfect in the regime of interest, and 
show that the effect of baryons on the power spectrum is not 
well approximated by a change in fl m h. These differences are 
important in order to be able to fit the 2dFGRS data with 
both parameters. The shape of the power spectrum is pri- 
marily dependent on Q m h, and only weakly dependent on h. 
We have therefore chosen to fit Q m h rather than f2 m . Sim- 
ilarly, the strength of baryon oscillations depend primarily 
on f2{,/fi m . We have therefore calculated model spectra for 
cosmologies chosen on a 40 x 40 x 40 grid in (Q m h, f2f,/fi m , h) 
space, covering 



0.1 < Q m h < 0.8, 
0.0 < Q b /Q m < 0.5, 
0.4 < h < 0.9. 



(•>) 



These limits effectively act as uniform prior probability den- 
sities for the parameters. A scale-invariant n — 1 primordial 
spectrum was assumed initially, following results from CMB 
analysis (e.g. Jaffe et al. 2000). We have numerically con- 
volved these power spectra on the scales 0.02 /iMpc -1 < 
k < 0.15 foMpc -1 usi ng th e window function fitting formula 
presented in Section 3.2. The normalization of the model 



spectra was allowed to vary to account for an unknown large- 
scale linear bias. 

Following analysis of the data using these model spec- 
tra, we have also created models on a finer 40 x 40 x 40 
grid covering 0.1 < Q, m h < 0.3, 0.0 < Qb/&m < 0.4 and 
0.4 < h < 0.9 in order to further constrain the fit in this 
region of parameter space. 



We emphasise that f2 m is the total matter density pa- 
rameter, i.e. Q m = Qcdm + Qb- Our results are not strongly 
dependent on tt v . 



6.2 Results 

The likelihood of each model has been estimated using a 
covariance matrix calculated from Gaussian realizations of 
linear density fields as in Section ^. Results presented in 
this Section were calculated using the covariance matrix es- 
timated by the method described in Section ^ for which 
Xmin = 27.5, given an expected value of 28. The best fit 
power spectrum parameters are only weakly dependent on 
this choice: we have tested a number of covariance matrices 
calculated for different input power spectra and find Xmin 
close to the expected value for all of them. In fact, we used 
an iterative procedure leading to this choice of covariance 
matrix. Initially we adopted a Q m h = 0.25, f2f,/f2 m = 
power spectrum, and then adopted the values Q m h = 0.2 
and Qb/Qm = 0.16 which are close to the best fit values 
determined with this covariance matrix. 

The likelihood contours in Qb/Q m versus Q m h for this 
fit are shown in Fig. ^. At each point in this surface we have 
marginalized by integrating the likelihood surface over the 
two free parameters, h and the power spectrum amplitude. 
The result is not significantly altered if instead, the modal, 
or Maximum Likelihood points in the plane corresponding 
to power spectrum amplitude and h were chosen. The like- 
lihood function is also dependent on the covariance matrix 
(which should be allowed to vary with cosmology) , although 
the consistency of result from covariance matrices calculated 
for different cosmologies shows that this dependence is neg- 
ligibly small. Thus C oc exp(— X 2 /2) in practice. 

Fig. [5] shows that there is a degeneracy between Q, m h 
and the baryonic fraction Qb/Sl m . However, there are two 
local maxima in the likelihood, one with £l m h ~ 0.2 and 
~ 20% baryons, plus a secondary solution Q m h ~ 0.6 and 
~ 40% baryons. 

Assuming a uniform prior for h over a factor of 2 is 
arguably over-cautious, and we have therefore multiplied by 
a Gaussian prior h = 0.7±10% in Fig. [| This corresponds to 
multiplying by the likelihood from external constraints such 
as the HST key project (Freedman et al. 2000). The effect is 
to tighten the contours around the above two models. The 
low-density model now becomes approximately 



n m h = 0.20 ± 0.03; n b /n m = 0.15 ±0.07. 



(6) 



The errors quoted are rms errors, and have been calcu- 
lated by integrating over the branch of solutions of interest. 
Analysing mock catalogues drawn from the ACDM Hubble 
volume simulation produces similar rms errors, and shows 
that, for each parameter, the interval defined by the ap- 
propria te e rror is close to a 68% confidence interval (see 
Section 



6.4) 



The 2dFGRS data are compared to the best-fit linear 
power spectra convolved with the window function in Fig. ^. 
This shows where the two branches of solutions come from: 
the low-density model fits the overall shape of the spectrum 
with relatively small 'wiggles', while the solution at Q. m h ~ 
0.6 provides a better fit to the bump at k ~ 0.065 ftMpc -1 , 
but fits the overall shape less well. 
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Figure 5. Likelihood surfaces for the best fit linear power spec- 
trum over the region 0.02 /iMpc" 1 < k < 0.15 ft, Mpc -1 . The 
normalization is a free parameter to account for the unknown 
large-scale biasing. Contours are plotted at the usual positions for 
one-parameter confidence of 68%, and two-parameter confidence 
of 68%, 95% and 99% (i.e. -21n(£/£max) = 1,2.3,6.0,9.2). We 
have marginalized over the missing free parameters (h and the 
power spectrum amplitude) by integrating under the likelihood 
surface. 



LO 
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"m h 

Figure 6. Likelihood surfaces for the best fit linear power spec- 
trum over the region 0.02 h Mpc -1 < k < 0.15 /iMpc -1 , as in 
Fig. ^, but now adding a prior on h: h = 0.7 ± 10%. This tight- 
ens the constraints. This result is compared to estimates from 
X-ray cluster analysis (Evrard 1997), big-bang nucleosynthesis 
(O'Meara et al. 2001) and recent CMB results (Netterfield ct al. 
2001; Pryke et al. 2001). Note that we have plotted the CMB 
result following the reasonable approximation that f!j,/i 2 and 
Sl ct l m /i 2 were independently determined by each of these anal- 
yses. 



k ranj 


;e / ft, Mpc -1 


n 


assumed 


ML parameters 


min 


max 




!!m for r(z) 


flmh 




0.02 


0.15 


1.0 


1.0 


0.23±0.03 


0.18±0.07 


0.02 


0.15 


1.0 


0.4 


0.20±0.03 


0.16±0.07 


0.02 


0.15 


0.9 


0.3 


0.22±0.03 


0.12±0.07 


0.02 


0.15 


1.1 


0.3 


0.18±0.03 


0.19±0.07 


0.015 


0.15 


1.0 


0.3 


0.20±0.03 


0.14±0.07 


0.03 


0.15 


1.0 


0.3 


0.20±0.03 


0.15±0.07 


0.02 


0.10 


1.0 


0.3 


0.17±0.04 


0.18±0.08 


0.02 


0.12 


1.0 


0.3 


0.18±0.03 


0.17±0.07 


0.02 


0.15 


1.0 


0.3 


0.20±0.03 


0.15±0.07 



Table 1. Maximum Likelihood Q m h and f2f,/f2 m parameters for 
fits to the 2dFGRS power spectrum, varying the range of fc-space 
fitted, the power-law index n of the primordial spectrum, and the 
matter density of the flat cosmology assumed to estimate the co- 
moving distance to each galaxy. All of these fits used a covariance 
matrix calculated from Gaussian realizations of a H m /i = 0.2, 
f2(,/f2 m = 0.16 CDM power spectrum. 




6.3 Robustness of the fit 

We have tried varying the range of k for the fit, the assumed 
power-law index of the primordial fluctuations, and the as- 
sumed geometry. The best fit Sl m h and Q.^/Q. m for a variety 
of assumptions are presented in Table |l| along with approx- 
imate errors, h = 0.7 ± 10% was assumed for this analysis. 
The shape of the likelihood surfaces and the position of the 
minimum recovered from each of these fits are similar, and 
the ML values generally change by <C la. The main effect of 
changing assumptions is in how rapidly the likelihood falls 
away from the ML point. 

Perhaps the main point to emphasize here is that the 



Figure 7. The 2dFGRS data compared with the two preferred 
models from the Maximum Likelihood fits convolved with the 
window function (solid lines). Error bars show the diagonal ele- 
ments of the covariance matrix, for the fitted data that lie be- 
tween the dotted vertical lines. The unconvolved models are also 
shown (dashed lines). The U m h ~ 0.6, fl b /U m = 0.42, h = 0.7 
model has the higher bump at k ~ 0.05 h Mpc -1 . The smoother 
Q m h ~ 0.20, fl b /Q m = 0.15, h = 0.7 model is a good fit to the 
data because of the overall shape. 
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results are not greatly sensitive to the assumed tilt of the 
primordial spectrum. We have used the CMB results to mo- 
tivate the choice of n — 1, but it is clear that very substan- 
tial tilts are required to alter our conclusions significantly: 
n ~ 0.8 would be required for the zero baryon model to 
become an acceptable fit, within la of the preferred model. 



6.4 Fitting to mock data 

Using the technique described in Section we have tried 
to recover the input parameters of the ACDM simulation, 
Q m h — 0.21 and fi;,/fi m — 0.13, from 35 redshift-space 
galaxy catalogues drawn from this simulation to mimic the 
2dFGRS data. Having calculated power spectra for the cat- 
alogues (presented in panel c of Fig. we have fitted the 
data with models using a covariance matrix calculated as in 
Section ^[ The Maximum Likelihood (and secondary max- 
ima if present) fl m h and Qt/Qm parameters recovered from 
these catalogues are presented in Fig. [| There is a degener- 
acy between recovered parameters: the data trace a pattern 
similar to that of power spectra with similar shape in the re- 
gion 0.02 /iMpc -1 < k < 0.15 ftMpc" 1 . This pattern is sim- 
ilar to that determined from the 2dFGRS data (Fig. ^|). The 
degeneracy between Q. m h and Q,b/£l m is weakly broken with 
the models favouring approximately the correct baryon frac- 
tion. For each maximum within the region 0.1 < Q. m h < 0.3 
& 0.0 < fii,/fi m < 0.4, we have integrated over the likeli- 
hood and find rms values similar to those obtained from the 
2dFGRS data. Of the 35 catalogues modelled, 27 (77%) have 
likelihood maxima within 1 rms of the true Q m h value, (i.e. 
they have likelihood maxima with 0.18 < Q, m h < 0.24). We 
also find that 24 (69%) have likelihood maxima within 1 rms 
of the true baryon fraction (i.e. they have likelihood maxima 
with 0.06 < flb/Qm < 0.20). This gives us confidence that 
the quoted errors for the best-fit parameters derived from 
the 2dFGRS data are sound. 

Combining the likelihood surfaces calculated from the 
13 non-overlapping catalogues results in best-fit parame- 
ters Q m h = 0.20 ± 0.15 and VL b /Q, m = 0.10 ± 0.05, with 
— 2 ln(£true/£ max ) < 1.0 for the combined likelihood. This 
offset is statistically acceptable, and suggests that, if there 
is a systematic bias in determining Q m h and fij,/n m from 
these catalogues, it is at a level well below the errors on the 
recovered parameters from any single catalogue. 



7 CONCLUSIONS 

We have shown that the present 2dFGRS data allow the 
galaxy power spectrum to be measured to high accu- 
racy (10-15% rms) over about a decade in scale at k < 
0.15 /iMpc -1 . We have carried out a range of tests for sys- 
tematics in the analysis and a detailed comparison with re- 
alistic mock samples. As a result, we are confident that the 
2dFGRS result can be interpreted as giving the shape of the 
linear-theory matter power spectrum on these large scales, 
and that the statistical errors and covariances between the 
data points are known. 

By fitting our results to the space of CDM models, we 
have been able to reach a number of interesting conclusions 
regarding the matter content of the universe: 
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Figure 8. 2D likelihood maxima in the (Q m h, Q(,/f2 m ) plane cal- 
culated from 35 redshift-space catalogues drawn from the ACDM 
Hubble volume simulation. We have marginalized (by integrat- 
ing under the likelihood) over h and P(k) amplitude assum- 
ing a Gaussian prior h = 0.7 ± 0.07 and a uniform prior for 
the amplitude. Primary maxima are given by the solid circles, 
other maxima with — 2 ln(£/£ max ) < 9.2 (of which there were 
~ 1.7 for each catalogue) by the open circles. Primary max- 
ima for which there was no secondary maximum are shown by 
the large solid circles. The locus of recovered parameters has a 
similar pattern to that of power spectra with similar shape over 
0.02/iMpc -1 < k < 0.15 hMpc -1 , and is also similar to that 
seen in the likelihood contours deduced from the 2dFGRS data. 
Solid lines mark the input parameters for the simulation. 



(1) The power spectrum is close in shape to that of a Q. m h — 
0.2 model, to a tolerance of about 20%. 

(2) Nevertheless, there is sufficient structure in the P(k) 
data that the degeneracy between Qb/Sl m and Q m h is 
weakly broken. The two local likelihood maxima have 
(tt m h,n b /fl m ) ~ (0.2,0.15) and (0.6,0.4) respectively. 

(3) The evidence for detection of baryon oscillations in 
the power spectrum is presently modest, with a like- 
lihood ratio of approximately 2.7 between the favoured 
model and the best zero-baryon model. Conversely, a 
large baryon fraction can be very strongly excluded: 
0.b/Q.m < 0.28 at 95% confidence, provided Q m h < 0.4. 

(4) These conclusions do not depend strongly on the value 
of h, although they do depend on the tilt of the primor- 
dial spectrum, with n ~ 0.8 being required to make a 
zero-baryon model an acceptable fit to the data. 

(5) The sensitivity to tilt emphasizes that the baryon sig- 
nal comes in good part from the overall shape of the 
spectrum. Although the eye is struck by a single sharp 
'spike' at k ~ 0.065 h Mpc -1 , the correlated nature of 
the errors in the P(k) estimate means that such fea- 
tures tend not to be significant in isolation. We note 
that the convolving effects of the window would require 
a very substantial spike in the true power in order to 
match our data exactly. Such over-fitting is not pos- 
sible within the compass of conventional models, and 
the conservative conclusion is that the apparent spike 
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is probably enhanced by correlated noise. A proper sta- 
tistical treatment is essential in such cases. 

It is interesting to compare these conclusions with other con- 
straints. Averaging the results of Netterfield et al. 2001 and 
Pryke et al. 2001, the current CMB data require Q m h 2 — 
0.15 ± 0.03, fl b h 2 = 0.0215 ± 0.0025, together with a power- 
spectrum index of n = 0.97 ± 0.06, on the assumption of 
pure scalar fluctuations. If we take h = 0.7 ± 10%, this gives 



Q m h = 0.21 ± 0.05; n b /Q, m = 0.14 ± 0.03, 



(7) 



in remarkably good agreement with the estimate from the 
2dFGRS 



Q m h = 0.20 ±0.03; 



0.15 ±0.07. 



(8) 



Latest estimates of the Deuterium to Hydrogen ratio in QSO 
spectra combined with big-bang nucleosynthesis theory pre- 
dict Q, b h 2 = 0.0205±0.0018 (O'Meara et al. 2001), in agree- 
ment with the latest CMB results. The confidence interval 
estimated from the 2dFGRS power spectrum overlaps both 
regions. X-ray cluster analysis predicts a baryon fraction 
O b /Q m = 0.127 ± 0.017 (Evrard 1997) which is again within 
la of our preferred value. 

The above limits are all shown on Fig. ^| and paint a 
picture of impressive consistency: it appears that we live 
in a universe that has Q m ~ 0.3 with a baryon fraction 
of approximately 15%. The precision of this statement will 
improve greatly with completion of the 2dFGRS. Doubling 
the sample size will improve the errors on the baryon fraction 
by much more than ~ V2- The window function will be more 
compact, so the signatures of baryon oscillations should be 
seen very clearly - or we will see that some non-standard 
alternative is required. 
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